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We study the existence of self-dual nontopological vortices in generalized Maxwell-Higgs models 
recently introduced in Ref. [3]. Our investigation is explicitly illustrated by choosing a sixth-order 
self-interaction potential, which is the simplest one allowing the existence of nontopological struc¬ 
tures. We specify some Maxwell-Higgs models yielding BPS nontopological vortices having energy 
proportional to the magnetic flux, >l>s, and whose profiles are numerically achieved. Particularly, 
we investigate the way the new solutions approach the boundary values, from which we verify their 
nontopological behavior. Finally, we depict the profiles numerically found, highlighting the main 
features they present. 


I. INTRODUCTION 

This work deals with the presence of nontopological 
vortices in generalized Maxwell-Higgs models recently 
introduced in Ref. [I|, whose dynamics is controlled 
by two positive functions depending on the Higgs field 
only. Before focusing on this, however, it seems inter¬ 
esting to first comment on topological solutions. As 
one knows, topologically nontrivial configurations, gener- 
ically named topological defects, are frequently described 
as static solutions to some nonlinear classical field models 
allowiM for the spontaneous symmetry breaking mecha¬ 
nism (2|. In this sense, such structures have been studied 
intensively all over the years, their consequences being 
applied to many areas of interest, specially in Cosmol¬ 
ogy, since they are known to be formed during symmetry 
breaking phase transitions Q. 

The interesting point to be raised is that, according 
to the usual approach, these solutions can be found by 
solving a given set of first-order differential equations 
(instead of the second-order Euler-Lagrange ones) 0. 
The main advantage of following such a prescription is 
that it also reveals that the resulting solutions are sta¬ 
ble against decaying into the respective mesons (i.e. they 
saturate a lower bound for the energy of the model), with 
their energy depending only on the physical fields bound¬ 
ary conditions (i.e. the topology). The simplest solu¬ 
tion is the static kink coming from a (l-l-l)-dimensional 
model described by one real scalar field only @ . In addi¬ 
tion, there are also vortices and magnetic monopoles: the 
first ones are born in (l-|-2)-dimensional gauged Abelian- 
Higgs theory while the second ones stand for the 
static solutions appearing in (l-|-3)-dimensional gauged 
non-Abelian-Higgs models Q- 

Recently, many authors have also investigated topolog¬ 
ical configurations coming from generalized models pos¬ 
sessing nonusual dynamics. In these works interesting 
new properties have been discovered, including variations 
on the shape of the resulting solitons, see Refs. [ 1 , 0 . 
In addition, some of these models have been verified to 
support first-order solutions and a lower bound for their 


energy, see Ref. 0- Here, it is important to point out 
that the main motivation regarding noncanonical dynam¬ 
ics comes in a rather natural way from the string theory, 
and that interesting applications of such an idea have 
been found in connection with the accelerated inflation¬ 
ary phase of the univers e flTI . tachyon matter [T^, dark 
matter [T^, and others |l4|. 

In field theory models, finite energy solutions require 
that the self-interacting potential goes to zero as the 
field solutions approach their asymptotic profiles. Stable 
nontopological solutions satisfy asymptotical boundary 
conditions that imply finite energy and null topological 
charge (becoming stable due to conservation of Noether 
charges). Nontopological defects with vorticity n > 0 
have null profiles at origin and far from it, requiring 
a 4>^ potential, as it occurs in the Chern-Simons-Higgs 
model [l^ and in some extended electrodynamics (iq . 
Nontopological solutions were initially studied in models 
constituted of one and two scalar fields, with applications 
in QCD see also Ref. [l^. Some investigations 

about the conditions under which modified Lagrangians, 
composed of generalized kinetic terms, yield nontopolog¬ 
ical soliton solutions were undertaken in Ref. [^ . 

The impossibility to attain nontopological solutions in 
a Maxwell-Higgs model is associated to the usual self¬ 
interacting potential, which provides BPS solutions. A 
generalized Maxwell-Higgs model with dielectric func¬ 
tions depending on the scalar field, was recently pro¬ 
posed [I| , and may be used to achieve nontopological so¬ 
lutions. In the present work, we investigate new aspects 
of the Maxwell-Higgs model of Ref. [l|, showing that 
it also supports nontopological configurations whenever 
supplemented with self-interacting potentials that engen¬ 
der symmetric minimum. 

In order to present our results, this work is orga¬ 
nized as follows. In Sec. H, we review the way the 
adopted nonusual Abelian-Higgs model engenders self¬ 
duality, embracing the first-order equations and the lower 
bound for the overall energy (defined in terms of an aux¬ 
iliary function conveniently introduced). In Sec. HI, we 
define the scenarios we consider by choosing the Higgs 
potential suitably. In addition, we characterize our so- 
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lutions by investigating the way they approximate the 
boundary values, which the physical fields are supposed 
to reach (near the origin and asymptotically). We also 
show that, in the appropriate limit, these models ad¬ 
mit the very same analytical solutions and related con¬ 
clusions. We go further by presenting numerical pro¬ 
files, which explicitly describe the general properties of 
the solutions. We still point out the existence of models 
engendering entirely analytical nontopological solutions 
(i.e. not approximate ones). Finally, in Section IV, we 
present our concluding remarks and perspectives regard¬ 
ing future investigations. 


II. THE MODEL 


We begin by reviewing the first-order framework firstly 
introduced in Ref. [I| for attaining topological solutions, 
whose starting-point is the (I -I- 2)-dimensional general¬ 
ized Maxwell-Higgs Lagrangian density, 

^ + w (|<(.|) - Um , (1) 

where = d^A,j — is the electromagnetic field 

strength tensor, is the minimal co¬ 

variant derivative. Here, we use standard conventions, in¬ 
cluding the plus-minus signature for the Minkowski met¬ 
ric. The positive functions, G'(|())|) and w{\(j)\), general¬ 
ize the dynamics of the usual Maxwell-Higgs model. It 
is easy to show that the Gauss law of this model is satu¬ 
rated by the temporal gauge, Hq = 0, so it describes only 
magnetic configurations, as the usual case. 

Our aim is looking for radially symmetric configura¬ 
tions coming from ©• For such purpose we implement 
the standard stationary vortex Ansatz, 


(r, 6) = vg (r) e*”®, A (r, 9) =-—{a{r) - n), (2) 

er 

where n = ±1,±2,±3... stands for the winding number 
of the topological configuration. 

We focus our attention on the self-dual configura¬ 
tions of Lagrangian o, governed by two coupled first- 
order equations that minimize the total energy of the 
model. In order to obtain such equations, we follow the 
usual Bogomornyi-Prasad-Sommerfield (BPS) formalism 
by writing the energy density of this model in the radially 
symmetric Ansatz, 


e 




/1 da 
\r dr 


2 


-I- v'^w 



(3) 


where U = U (g). In the present scenario, self-duality 
only holds when G(|(/)|), rc(|(/)|) and U {\(j)\) are con¬ 
strained to each other by 

— V GU = —\/2ev^wg. (4) 

dg 


Under such a constraint, the energy density Eq. m can 
be rewritten as, 


e = 



+ VW 


da ag 
dr r 


- — — (aV2Gu) 
er dr \ J 


(5) 


from which one concludes that the corresponding total 
energy saturates its lower bound when g (r) and a (r) 
satisfy 


dg 

dr 

_ ,ag 

r 

(6) 

1 da 


(7) 

er dr 

= 


These are the BPS or self-dual first order equations of the 
generalized model. Saturating the first-order equations, 
the energy density is shown to be 

1 dH 

£bps — T 1 1 (8) 

er dr 

where we have defined the auxiliary function 

H (r) = ay/2GU. (9) 


By using BPS equations, we also can write the correspon¬ 
dent energy density as 


ebps = 2U + 2v^w(^'^y , ( 10 ) 

which will be positive definite for w{g) > 0. In summary, 
given G, w and U constrained by relation dH), the first- 
order equations ([61) and © give rise to self-dual solutions 
for both g(r) and a(r) which satisfying the boundary 
conditions. 


5(0) = 0, g{oo) = 1, 

(11) 

a (0) = n, a (oo) = 0, 

(12) 


engendering topological configurations. These solutions 
possess finite energy obtained by integrating Eq. ([5]) over 
the plane, i.e. 

Ftps = ^TT y rSbpsdr = ±2tt [H (0) - H (oo)], (13) 

with H (0) and H (oo) = 0 being the values of F[{r) at 
the boundaries. In this model, suitable choices of func¬ 
tions G and w allowed to achieve finite energy BPS so¬ 
lutions associated with potentials U different from the 
usual (jA one [l|, with the energy resulting proportional 
to the magnetic flux. So far, no study about the possi¬ 
bility of attaining BPS nontopological solutions in this 
Maxwell-Higgs framework was developed. 
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III. NONTOPOLOGICAL VORTICES following behavior at boundaries, i.e. when g (r) —>■ 0, 


In this Section, we show that the generalized Maxwell- 
Higgs model of Ref. [l| may also support nontopologi- 
cal radially symmetric solitons, assuming some specific 
choices for the generalized functions G and w, and that 
the first-order equations dH]) and © are solved fulfilling 
suitable boundary conditions. We can show that the en¬ 
ergy of these BPS nontopological magnetic vortices is 
proportional to the magnetic flux not necessarily 
quantized in this case. From now on and for simplic¬ 
ity, we choose e = v = 1, and consider only the case 
n > 0, corresponding to the upper signs in Eqs. ®, 0 , 

®, CSl). 

In order to obtain nontopological configurations, the 
profile functions g (r) and a (r) must obey the following 
boundary conditions for n ^ 0 , 


II 

0 

g 

II 

0 

(14) 

a ( 0 ) = n, a(oo) = —an, 

(15) 


where an stands for a positive real number calculated 
numerically. Under such boundary conditions, the mag¬ 
netic flux <I>B is immediately obtained 


<^( 5 ) = ^+70 + (18) 

9 

with the constants 7 , 70 , 72 , •.• characterizing such be¬ 
havior. The case 7 7 ^ 0 allows to express the total BPS 
energy as being proportional to the magnetic flux. In the 
following we fix 7 = 1 . 

By using the Eqs. (II3 and (1181) . we should investigate 
the way the Helds g (r) and a (r) behave near the bound¬ 
aries. We thus solve the first-order equations ® and © 
around the boundary values (ITT)) and (fT51) . obtaining the 
behavior of profiles near the origin. 


g(r) s 

- Cor", 

(19) 

a(r) s 

.. „ ^'-"0 2n+2 

2 (n + l) 

( 20 ) 


and asymptotically, 


9{i") 

a{r) 



— Ctn + 


1 

2 {an — 1 ) 


( 21 ) 

( 22 ) 


= 27r j rB (r) dr = 2 tt (n + an), (16) 

and is not necessarily quantized (i.e. an is not an inte¬ 
ger). 

Notwithstanding the magnetic flux be given by (Ull), it 
is still possible to define classical field models for which 
the magnetic flux defines the lower bound for the BPS 
energy (HSl). Note that this occurs only when the function 
H satisfies boundary conditions proportional to the ones 
of Eq. (I15L which requires a particular choice for the 
functions G and w. 

It is well-known that nontopological solutions arise in 
field theories, which possess at least one symmetric vac¬ 
uum, and asymptotical conditions compatible with finite 
energy. In canonical field models, the simplest Higgs- 
potential consistent with this requirement is the usual 
sixth-order one, 


U{g) = ^g^{l-g^)\ (17) 


Here, Cq and Gao are positive real constants obtained 
numerically by requiring proper behavior near and far 
from the origin, respectively. 

The point to be raised here is that the solutions (IT^ 
and o stand for a typical behavior of a topological 
vortex near the origin whereas Eqs. (1211) and (I22|) encode 
the typical nontopological profile far from the origin. It 
is worthwhile to say that in Ref. di , it was shown that 
the very same dependence in r (despite some numerical 
factors) characterizes the nontopological vortices of the 
usual Chern-Simons-Higgs model. 

For small g, replacing expressions 013 and ( 11 ® in the 
self-dual Eqs. © and ©, one achieves a Liouville’s equa¬ 
tion for g (r) (in accordance with the procedure of Ref. 
d^), whose analytical solution is 


9{r) 


2 (n + l) 
roVX 




2(n+l) 


+ 1 


(23) 


where A stands for the coupling constant for the scalar- 
matter self-interaction (supposed to be dimensionless and 
positive). In this case, the symmetric vacuum is g = 0, 
the asymmetric one standing for g = 1. Here, we point 
out that the simplest self-dual Maxwell-Higgs model (sat¬ 
urated by the usual | 0 |'^-potential) presents only the 
asymmetric vacuum g = 1 and does not support nontopo¬ 
logical structures. Within our generalized scenario, we 
circumvent this problem by adopting the |(/)|®-potential 
(II3, whilst choosing the functions G and w conveniently. 

Without loss of generality, and in the context of La- 
grangian ©. we suppose that the function G {g) has the 


with its maximum located at r = ro (n/ (n -|- 2 ))^^^^"'“''^^ 
(for n sufficiently large, its location matches tq). The cor¬ 
responding expression for the gauge function a (r) reads 
as 


a (r) = —n — 2 + 


2 (n + l) 



2(n-|-l) 


+ 1 


(24) 


The (approximate) solutions attained this way fulfill con¬ 
ditions (fTTl) and (TTSl) . and hold for r —>■ 0 or r —>■ 00 , with 
an = n + 2, confirming the nontopological behavior. 
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IV. SOME MODELS ENDOWED WITH 
NONTOPOLOGICAL SOLUTIONS 

We now show some models belonging to the generalized 
Maxwell-Higgs Lagrangian o, endowed with the sixth- 
order potential da- These models are better character¬ 
ized by choosing a specific function w ( 5 ), from which the 
constraint dH) yields G {g) behaving as (fT51) . 

The first model is defined by the potential (HZD and the 
functions 


( 5 ) = ( 5 ^-b 1) , G (g) = . (25) 

From Eqs. and 0 , we get the self-dual equations of 
this model, 

^ ^ ^ ^ 3Ag^ (g^ - l) 

dr r ^ r dr g^ + 3 ’ 

which was firstly studied in Ref. [l[ (except for some nu¬ 
merical factors) for analysis of topological solitons. The 
BPS nontopological configurations of this specific model 
satisfy the boundary conditions (fldl) and CSl), with en¬ 
ergy given by Eq. (US. Erom the choices of this model, 
we calculate H (0) = nX and H ( 00 ) = —anX, so that the 
total energy is 


At this point, it is interesting to note that because of 
Eqs. (E51) and (1^ the values of the winding number n 
are limited by the requirement g < 1 , otherwise we do 
not obtain solutions satisfying the boundary conditions 
or with finite energy. 

The third |(()|®-model to be scrutinized is one that af¬ 
fords equations, which have an analytical solution. It is 
specified by the functions 

w = 2X{l-g^), G{g)= ~ , (30) 
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providing the following self-dual equations: 


dg ag 1 da 
dr r ’ r dr 


-Xg\ 


(31) 


The Eqs. (iMl) are exactly solvable, having analytical 
solutions given by (|23)l and (l24)l . which hold for all r, 
not only at origin or at infinity. These solutions fulfill 
conditions da and la, with an = n + 2, possessing 
magnetic flux <i)s = 47 r (n -b 1 ) and total energy Ftps = 
47rA (n -b 1 ). 

The requirement of a positive-definite energy density 
(fTOl) or positive w (g) in (l30l) is equivalent to demanding 
g < 1 , which holds only when 


Xrl > (n -b (32) 


Ftps = 2 ttX (n + an) ■ (27) 

thus we verify that the energy of the resulting self-dual 
structures is proportional to their magnetic flux. 

The second |(/)|®-model we study is defined by the func¬ 
tions 


w{g) = X, G{g) = -^, (28) 
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whilst the self-dual equations are 


dr 


7 

r 


1 da 
r dr 


Xg^ (g 2 - 1) , 


(29) 


whose topological solutions were analyzed in Ref. [l|. 
We observe that Eqs. (l29l) mimic those of the standard 
self-dual Chern-Simons-Higgs model m for A = 2 /k^ (k 
standing for the Chern-Simons constant). Despite satis¬ 
fying the very same differential equations, the solutions 
we describe in this work are physically different, since 
they present no electric charge. 

By solving the BPS equations (l29ll around the bound¬ 
ary values (HI and m, we verify that g (r) and a (r) 
are given as in Eqs. (HU and (I20p near the origin, while 
their asymptotic behavior coincides with the ones of (l2T]l 
and (1^ . Thus, the solutions coming from (1^ indeed 
stand for nontopological self-dual vortices possessing no 
electric charge. The total BPS energy is computed by us¬ 
ing Eqs. (fT71) . (051) . (fTTl) and (ITSl) . leading to H (0) = nX, 
H ( 00 ) = —a„A, and F^ps = 27rA (n + a„) = A4 >b, being 
proportional to the magnetic flux again. 


is satisfied. Therefore, for a given value of the coupling 
constant A, the possible values of the vorticity are limited 
in such a way that larger values of rg imply larger values 
of n. 

We now show the numerical solutions for fixed A = 100 
and ro = 1. For these values, the winding number is 
restricted to 1 < n < 4, 1 < n < 5, or 1 < n < 8 , for 
the first, the second and the third models, respectively. 
It is worthwhile to note that a similar situation holds 
in the Chern-Simons-Higgs model [l^, whose the BPS 
equations are the very same as those of the second model, 
see (1551) . 

In Figs. 1 and 2, we depict the solutions for g (r) 
and a(r), respectively. The first ones are peaked at 
rings around the origin, whereas the second ones stand 
for peaks centered at the origin. It is interesting to 
note that, for the first model one finds ai = 3.03795, 
a 2 = 4.07752, and ai = 3.02977, 02 = 4.06176, for the 
second model. Therefore, we find a small deviation of 
the value an = n -b 2 . 

The Figs. 3 and 4 show the solutions for the magnetic 
field B {r) = —{da/dr)/r and for the self-dual energy 
density (O, respectively. In general, all the profiles ex¬ 
hibit the same behavior, being ring-like peaks centered 
at the origin, and vanishing asymptotically. 


V. ENDING COMMENTS 

In this work, we have investigated the way the gen¬ 
eralized Maxwell-Higgs model o gives rise to self-dual 
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r 


FIG. 1: Solutions to g {r) calculated via (i) (1261) (solid blue 
line for n = 1, dotted blue line for n = 2), (ii) (1291) (long- 
dashed red line for n — 1, dash-dotted red line for n = 2), 
and (iii) (EIJ (i.e., the analytical profiles (1231) and (1241) ; dashed 
black line for n = 1, space-dotted black line for n — 2). We 
have used A = 100 and ro = 1. 




r 


FIG. 4: Solutions to etps- Gonventions as in Fig. 1. 

nontopological magnetic vortices possessing finite energy 
and no electric charge. In order to present our results, 
we have first reviewed the basic features of the model, 
which engenders self-duality only when a particular con¬ 
straint dH) is satisfied. We have specified our Maxwell- 
Higgs models by choosing the usual sixth-order potential 
(II3 and changing the function w (|^|) conveniently, from 
which we have obtained the corresponding expression for 
G {\(j>\) from the constraint (|4|). Without loss of gener¬ 
ality, we have normalized the function G{g) in such way 
that the energy of the nontopological solutions is propor¬ 
tional to the corresponding magnetic flux. 

The BPS equations of the models here analyzed when 
g —?> 0 reduce to exactly the equations (IdTll . In this held 
regimen the solutions are exact and given by Eqs. (1251) 
and (1251) . These equations and their solutions are invari¬ 
ant under the following scaling transformation, 

r —>• (5r , g {6r) = 6~^g (r), a (Sr) = a{r), (33) 


FIG. 2: Solutions to a{r). Conventions as in Fig. 1. 



r 


FIG. 3: Solutions to B (r). Conventions as in Fig. 1. 


with 6 the transformation parameter. Such an invariance 
enhances the solutions of the models because there must 
exist a free parameter dehning a family of inhnite solu¬ 
tions. In |(/)|®-models dehned in (I -|- 2)-dimensions the 
A coupling constant is dimensionless. Thus it becomes a 
free parameter characterizing the solutions. For example, 
A determines the behavior of the solutions of m given 
by Eqs. (1251) and (125)) when g ^ 0. This way, when 
r —>• 0: 


gir) 


2 (n+l) 

roVX Vo/ 


and for r —>■ oo: 


9{r) 


2 (n -I-1) /ro \ "+2 

roVX + 


(34) 


(35) 


An important result is that for Hxed A and tq, we can 
determine the values of n for which exist nontopological 
BPS vortices, as given by Eq. (1521) . In fact, we have 
verified that the winding number values are limited by 
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the condition g <1. Despite this inequality (l32l) is exact 
only for the solution of the set of BPS equations (ISTl) and 
it gives an upper-bound for the values of n in all models 
we have considered. 

An other result is related to the BPS energy of the 
nontopological vortices; for the models we have studied, 
the energy depends explicitly on the free parameter A, 
see (l?7l) . 

The stability issue of the new nontopological BPS so¬ 
lutions is under investigation and the results will be re¬ 
ported elsewhere. 


An other interesting issue concerns the search for BPS 
solutions in generalized Born-Infeld-Higgs models. This 
issue is now under consideration, and we expect interest¬ 
ing results for a future contribution. 
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